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INVESTIGATION OF WO BLUFF SHAPES IN AXIAL FREE FLIGHT 
OVER A MACH NUMBER RANGE FROM 0.55 TO 2.15 
By Lucille C. Coltrane 

SUMMARY 


A fineness -ratio-2. 71 right circular cylinder and a fineness-ratio- 
k.-O circular cylinder with a rounded nose and a flared afterbody have 
been tested in free flight over a Mach number range of O.J5 to 2.15 and 

a Reynolds number range of 1 X 10° to 12 x 10°. Time histories, cross 
plots of force coefficients, rolling velocity, and longitudinal -force 
coefficient are presented for both cylinders. In addition, cross plots 
of moment coefficients and plots of the normal -force curve slope and 
the aerodynamic center are presented for the fineness-ratio-2.71 cylinder. 
The average aerodynamic center of the right circular cylinder moved rear- 
ward with decreasing speeds until at the subsonic Mach numbers it remained 
approximately constant and comparisons of the drag data of this test with 
wind tunnel and other free-f light data show good agreement. An appreciable 
decrease in drag was observed when the data of the present test of the 
rounded nose cylinder were compared with data of a right circular cylinder 
of a similar configuration. 


INTRODUCTION 


The dynamic stability characteristics of shapes having very high 
drag and very low lift are of interest to designers of various aerodyna- 
mic stores (see refs. 1 and 2) and shapes such as a reentry missile. One 
configuration which meets these requirements is the circular cylinder in 
axial flow. An investigation is being conducted by the Langley Pilotless 
Aircraft Research Division with the rocket -boosted free -flight-model 
technique to obtain stability data on bluff shapes in axial free flight. 

The present investigation was conducted to extend the data in ref- 
erence 5 to higher Mach numbers and to determine the effect of a change 
in nose shape on the aerodynamic characteristics of a bluff shape. For 
this purpose flight tests were made of two additional cylinders, a right 
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circular cylinder with a fineness ratio of 2.71 of the same design as 
the fineness-ratio-2.56 cylinder of reference 5 and a circular cylinder 
with a rounded nose and a flared afterbody of the same design except for 
nose shape as the fineness-ratio- 4.0 cylinder of reference 5. The 
results of these flight tests are presented in this paper. These tests 
covered a Mach number range from 0.55 to 2.15 and a Reynolds number 

range from 1 x 10^ to 12 x 10^ based on cylinder diameter. The free- 
flight tests were conducted at the Langley Pilotless Aircraft Research 
Station at Wallops Island, Va. 


SYMBOLS 


The data are presented relative to the body axis system and the 
positive directions of the force coefficients, moment coefficients, and 
angular velocities are shown in figure 1 . The various symbols used 
throughout the paper are defined as follows : 


a accelerometer reading, g units 



longitudinal -force coefficient, 



m 


pitching -moment coefficient, 


— (0 - t<p) 
qSd 


j 





2 V 


C 


N 


normal -force coefficient, 




dC N 

da 


yaw ing -moment 

coefficient, 

iz.(* 

qSd 

+ 0<p) 

lateral -force 

coefficient. 

a. 

w/s 



t,cg 

q. 

diameter of cylinder, ft 





NACA RM L58A16 


5 


g 

I 


acceleration due to gravity, ft/sec^ 

O 

moment of inertia, slug-ft 1 ^ 


X X> X Y’ I Z 
k 

l 

M 
m 

q. 

R 
r 
S 
t 
V 
W 
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N) 

aX 


p 

<p 


moments of inertia, slug-ft^ 


radius of gyration, ft 

m 


length of cylinder, ft 
Mach number 
mass, slug 

dynamic pressure, lb/sq ft 
Reynolds number based on cylinder diameter 
radius of cylinder nose, in. 
cross-sectional area of cylinder, sq ft 
time, sec 

free -stream velocity, ft/sec 
weight of cylinder, lb 

distance along cylinder from nose, ft, or, when used with a 
subscript, distance from center of gravity, positive 
forward, ft 

angle of pitch 

nonrolling damping constant, l/sec . 
damping constant due to roll, l/sec 


relative -density factor, 

air density, slugs/cu ft 
angle of roll, radians ’ 
angle of yaw, radians 


4m 

pSd 
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CUq 

basic oscillation frequency, radians/sec 

Ad 

component of total pitch frequency resulting directly from 
roll, radians/sec 

0 - ijr<p 

8, 0 — a 

pitching acceleration, g , radians/sec 2 

x n ,2 ~ Xn,l 

**11, eg 

, . . a n,l x n,2 “ a n,2 x n,l . , 

normal acceleration, — 2 — — - — , g units 

Xn,2 - ^,1 

*• • • 

i|r + 09 

a t 2 “ a t 1 / P 

yawing acceleration, g — - — , radians/sec 

x t,2 - ^,1 

a t,cg 

a t,l x t,2 ‘ a t, 2 * t,l 

transverse acceleration, — - - - — , g units 

x t,2 * x t,l 

Subscripts : 

ac 

aerodynamic center 

eg 

center of gravity 

l 

longitudinal 

n 

normal 

t 

transverse 

tot 

total 

1 

forward end of cylinder 

2 

rear end of cylinder 


A dot above a symbol indicates time rate of change of symbol, for 
example, 9 = d0 /dt . 


MODELS 

The physical characteristics of the models are presented in table I, 
and drawings of the models are shown in figure 2. Photographs of the 
models and model booster arrangements are presented in figure 3. 
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The fineness -ratio -2 . 71 model (designated model I) was a right cir- 
cular cylinder with a diameter of 8 inches and contained two small rockets 
which were mounted normal to the longitudinal axis and ahead of the center 
of gravity to give a yaw disturbance. The center of gravity of this 
model was located at 3! . 3 percent of the cylinder length behind the nose. 
The fineness -ratio -4.0 model (designated model II) was a circular cylinder 
with a diameter of 8.0 inches, a 13 ° flared afterbody, and a rounded nose 
and contained no small rockets. The afterbody and nose lengths are shown 
in figure 2 and the contour ordinates for the nose are shown in table II. 
The center of gravity of this model was located at 24.5 percent of the 
model length behind the nose. 

The cylindrical portions of the models were constructed of steel and 
were covered with a fiber-glass plastic shell. The forward portion of 
model I and the nose of model II were machined from solid steel. The 
afterbody flare of model II was constructed of l ami nated wood. 


INSTRUMENTATION 


Model instrumentation consisted of an NACA six-channel telemeter 
which transmitted data from six accelerometers located as follows : one 

normal and one transverse accelerometer in the forward end of the 
cylindrical portion of the model, one normal and one transverse accel- 
erometer in the rear of the cylindrical portion of the model, and two 
longitudinal accelerometers, one for high range and one for low range, 
behind the center of gravity. A measure of the signal strength trans- 
mitted from the loop antenna provided an indication of the roll rate of 
the cylinders since the strength of the signal varied with the cylinder 
roll position. 

Ground instrumentation included a CW Doppler radar unit to measure 
the velocity of the model, a modified SCR 5^4 tracking radar set to 
determine the flight path and a rollsonde receiver used as an additional 
measure of the rolling velocity. Fixed and tracking motion -picture 
cameras were used to observe the model during the first portion of the 
flight. Atmospheric data were obtained from a rawinsonde released immedi- 
ately before model flight. 


TESTS AND ANALYSIS 


Tests 


The models were ground launched at an angle of 70° from the 
zontal by means of a mobile launcher as shown in figure 3(c). A 
fuel Cajun rocket motor boos£gcl - the_m2j2£^^ to maximum velocity. 



hori- 

solid- 

Drag 
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flaps were incorporated, into the booster to increase the separating force 
between the model and booster at booster burnout. For model I, the drag 
flaps operated as designed, but, for model II, motion pictures showed the 
drag flaps opened prematurely, came off, and damaged the booster fins. 

The damaged fins probably increased the drag of the booster and, although 
separation was delayed about 0-3 second, a clean separation was experi- 
enced. In addition, the damaged booster fins caused the model-booster 
combination to experience a high rolling velocity immediately before sepa- 
ration. This high rolling velocity persisted throughout the flight of the 
model. Tracking radar showed that the models followed an approximately 
parabolic flight path. 


Analysis 

As was done in reference 4, the method of analysis used herein to 
determine the pitching-moment coefficients consists of using the data 
obtained from normal accelerometers located at two positions in the model. 
When the normal -force equation is combined with the equation describing 
the normal -accelerometer reading at any point along the x-axis of the 
model, it can be shown that 


g 


*n ,2 - Xn,l 


(%2 “ %l) = 9 - tf<P 


( 1 ) 


The pitching -mcment equation may be written as 


dSd^t = 1*0 - (IZ " I X )tf (2) 

if I xz = 0. Now, if Iy = Ig as is the case for the models, reported 
herein (see table I) , . . 


qSdC mtot = I y (0 - <pjr ) + I X <P^ 

and then, if Iy<p|r is small with respect to ly(0 - <pjr), 


Sot ■ 


(3) 


(4) 


For model I, equation (4) may be used directly since the roll rate 
was found to be small and the term I^cplr can neglected. 


For the aerodynamic -center location, where 


or 


d x cg 

ac dC N l l 


x ac 


dC n d + x cg 
dCy 1 Z 


the values for 
C m against 


dC m 


in 

and 

dC N 


and 

c n 


dC. 


dC ^- were taken directly from the polars of 
against Cy at Cjj or Cy near zero. 


In order to determine the normal -force coefficient curve slope as 
presented herein for model I and for the fineness-ratio-2.56 cylinder 
of reference 5, equation ( 42 ) in reference 5, which states that 




2 2 

^ + 


- AA 2 + 



( 5 ) 


was simplified by omitting the damping terms. Analysis of the data that 
were obtained after the firing of the two pulse rockets showed that the 
last three terms within the brackets were small when compared with the 
term. of the basic oscillation frequency cjq^. (See table III.) There- 
fore, they were neglected throughout the calculation. Then by the 
relation 


/ x 2 x l\ 

\ d d / dC m 

(k/d) 2 dC N 



5 

dC- 


( 6 ) 


it is possible to express 


in terms of 


as follows : 



The slopes S and — - were equal for the polars of model I. Since 

dCjj dCy 

variations of and ct were probably nonlinear, this average 

may be greater or less than values near Cjj = 0, depending upon the 
range of C N encountered in flight. 


An indication of the dynamic stability following the disturbance 
given the model by the pulse rockets was obtained through the use of 
equation (^ 4 - 3 ) in reference 5, where 



(8) 


ACCURACY 


The instrument inaccuracies, estimated to be ±2 percent of their 
calibrated range, are stated in coefficient form for representative Mach 
numbers as follows : 


Coefficient 

Fineness -ratio- 2 . 71 
model (model I) at 
Mach numbers of - 

Fineness -ratio - 4 . 00 
model (model II) at 
Mach numbers of - 

1.88 

1.26 

0.80 

2.00 

1.38 

0.80 

C N 

± 0.039 

± 0.094 

±0.250 

±0.038 

±0.085 

± 0.309 

c Y 

± 0.040 

± 0.095 

± 0.254 

± 0.039 

±0.087 

±0.316 

c c 

±0.051 

±0.030 

± 0.079 

±0.052 

±0.030 

± 0.104 
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PRESENTATION OF RESULTS 


The variation of test Reynolds number, based on cylinder diameter, 
with test Mach number is presented in figure 4. The model flight paths 
are presented as plots of altitude against horizontal distance in fig- 
ure 5 and the variation of velocity and dynamic pressure with t im e is 
shown in figure 6. 

Time histories of the normal -force coefficient, lateral-force coef- 
ficient, and Mach number are presented in figures 7 and 8 . Basic data 
cross plots of force and moment coefficients for model I are shown in 
figure 9* Examples of the variation of the normal -force coefficient with 
lateral-force coefficient for model II are shown in figure 10. An indi- 
cation of the rolling velocity against Mach number is shown in figure 11. 
For model I, as a function of Mach number, the variation of the average 
aerodynamic center in cylinder lengths behind the nose is presented in 
figure 12 , and the normal -force curve slope is presented in figure 13 . 

The measured longitudinal -force coefficient for both models is plotted 
against Mach number and is presented in figure 14. 


DISCUSSION OF RESULTS 


Time History 

A time history of M, C N , and Cy for model I is shown in fig- 
ure 7 . A disturbance at separation (t = 2.85) caused the instruments to 
exceed their calibrated ranges at each oscillation peak for about 1 sec- 
ond. Throughout the flight a coupled motion with respect to the body axis 
system was experienced by the model as the oscillations damped from large 
amplitudes at supersonic speeds to low amplitudes at subsonic speeds . 
Figure 8 presents a time history of M, Cjj, and Cy for model II. 


Basic Data Plots 

Plots of C N against Cy presented in figures 9 and 10 for various 

Mach numbers indicate the model motion. When the trim center is estimated, 
.the time history of the resultant force vector can be obtained and the 
; ,;arigular displacement of adjacent peaks gives an indication of the rolling 
'fstte. This use of the Cjj and Cy cross plots is discussed in detail in 

reference 5 . Also, as stated in reference 5 , whether the missile is above 
or below roll resonance is immediately determined, since the motion above 
resonance is characterized by inside loops, and the motion below roll 
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resonance by outside loops. Therefore, it can be seen that model I was 
below roll resonance and model II was above roll resonance; in fact, by 
comparing the plot of C N against C Y of figure 10 with a similar plot 

in figure 5 of reference 5, it is apparent that model II rolled at a rate 
of about three times the pitch frequency. This ratio is nearly the same 
as the ratio of the measured roll rate of model II and the nonrolling 
pitch frequency of a nearly similar fineness ratio 4.0 cylinder of ref- 
erence 3. Dashed lines in figure 9 indicate extrapolated values where 
the range of the instruments was exceeded but these data points were con- 
sidered accurate enough for obtaining an estimate of the roll rate. This 
estimate is shown with measured averages from the rollsonde in figure 11 
and good agreement is noted. 

Figure 9 shows that the variation of moment coefficient with force 
coefficient for model I was nonlinear in both the pitch and yaw planes, 
the greatest stable slope generally occurring at the zero force coeffi- 
cient. An attempt was made to determine the stability characteristics for 
model II from the disturbance received at separation; however, the scatter 
of data for the plots of C m against C N and C n against C Y pro- 
hibited the determination of a slope dC^JdC^ to be used in equation (6). 
Thus neither nor the aerodynamic -center location could be deter- 

mined. The relative importance of the various terms within the brackets 
of equation (5) for a model rolling about 100 radians/sec can be noted in 
table III for model II. 


Aerodynamic Center 

The variation of the aerodynamic center with Mach number for model I 
is shown in figure 12. The average aerodynamic center moved rearward as 
the speed decreased through the supersonic and transonic speeds but 
remained at a constant position of about 0.64 cylinder lengths behind 
the nose through the subsonic range. For comparison, a dashed-line curve 
faired from data of the fineness -ratio-2. 56 cylinder of reference 3 is 
also presented in figure 12. The comparison is considered good especially 
when the accuracies in the basic data are taken into account. 


Lift 

The variation of the normal -force curve slope with Mach number for 
model I is shown in figure 13- It is seen that Cjg decreased with 

CL 

decreasing speeds until at subsonic Mach numbers it remained approximately 
constant. Also shown in figure 13 is a curve of Cjj from faired data of 

the fineness -ratio-2. 56 model of reference 3 which shows good agreement. 


• • • 
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With known, was calculated to be -26 per radian at 

M = 1.03 and 22 per radian at M = O. 7 O. Since the term Aq is smal 1 

(see table III) and in these flight tests is difficult to obtain with 
good accuracy, these values should be considered as good estimates and 
not absolute values of the dynamic stability at these Mach numbers. 


Drag 

The variation of longitudinal -force coefficient with Mach number 
for the models of the present tests is presented in figure 14. The 
shaded areas on the figures indicate values of Cq obtained for all Cjj 

and Cy values of the tests. Also presented are the drag data for both 
models when Cjj < 0.1 and Cy < 0.1. For model I, comparison of these 

data with the drag of the fineness -ratio- 2. 56 cylinder of reference 3 
when Cjj < 0.1 and Cy < 0.1 and the drag data of reference 6 for the 

body of revolution where r = 0.00 inch shows good agreement. Compari- 
son of the data of model II with the data of the fineness -ratio -4.0 cylinder 
with flared afterbody of reference 3 when C N < 0.1 and Cy < 0.1 indi- 
cates that a considerable reduction in drag throughout the Mach number 
range, was obtained by the nose modification of the present test. The drag 
data of model II and the drag data of reference 6 (where r = 0,50 inch) 
show reasonable agreement at subsonic speeds . 


CONCLUDING REMARKS 


A fineness -ratio-2 . 71 right circular cylinder and a fineness -ratio - 
L.O cylinder with a rounded nose and a flared afterbody have been tested 
in rocket-boosted axial free flight over a Mach number range of O .35 to 

2.15 and a Reynolds number range of 1 x 10^ to 12 x 10^. When given a 
disturbance at supersonic speeds, the response of the f ine ness -ratio - 2. 71 
cylinder was a large -amplitude oscillation which damped to a low -amplitude 
sustained oscillation at subsonic speeds. The model experienced a low roll 
rate throughout the test. This value was less than the roll rate required 
for resonance. The moment -force curves were nonlinear, with greatest 
stable slope generally occurring at the zero force coefficient. The aver- 
age aerodynamic center moved rearward and the normal -force curve slope 
decreased with decreasing speeds until at the subsonic Mach numbers they 
remained approximately constant. Comparisons of the drag data of this 
test with wind-tunnel and other free -flight data show good agreement. 

The fineness -ratio -4.0 cylinder experienced a high rate of roll throughout 
the test. This value was considerably more than the roll rate required 
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for resonance. An appreciable decrease in drag was observed when the 
data of the present test of a rounded nose cylinder were compared with 
data of a right circular cylinder of a similar configuration. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., December 30, 1957* 
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TABLE I.- PHYSICAL CONSTANTS FOR CYLINDERS TESTED 


Constant 

Model I 

Model II 

W, lb 

97.25 

IO8.5 

Ij£, slug -ft ^ . . . 

O.I75 

0.238 

Iy, Slug-ft 2 . . . 

0.910 

1.430 

1%, slug-ft 2 . . . 

0.910 

1.430 

x C g/l ...... 

0.343 

0.245 

a, ft 

O.667 

O.667 

1, ft 

1.805 

2.667 

S, ft 2 

0.348 

0.348 


TABLE II.- NOSE ORDINATES FOR MODEL II 


0 

.01 
.10 
.20 
• 30 
.40 
.50 
.60 
.80 
1.00 

1.50 
2.00 

2.50 
3.00 
3.50 
4.00 


J, in. 


1.600 

2.198 

2,765 

2.964 

3.086 

3.177 

3.249 

3.309 

3.406 

3.483 

3.627 

3-733 

3.817 

3.887 

3.947 

4.000 
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TABLE III.- VALUES USED TO DETERMINE SK)PE OF 
NORMAL-FORCE COEFFICIENT 


Mach 

number 

^0 

A o 

AA 


em 

Model I 

1 

1.03 

22.5 

-1.33 

0 

2.75 

0.22k 

0.70 

i4.o 

+0.57 

0 

1.95 

0.22k 

Model II 

2.07 

37.0 

-2.33 

'1.13 

95*3 

0.190 
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(a) Model I. L- 57 - 1733 



(b) Model II. L-96363.I 

Figure 3*- Photographs of cylinders tested. 




(c) Model II on booster in launching position. 


L- 57-1356.1 


Figure 3.- Concluded. 
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0 .4 .8 1.2 1.6 2.0 2.4 


M 


Figure 4.- Reynolds number of tests based on cylinder diameter. 
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Horizontal distance, ft 


Figure 5- - 


Flight paths of cylinders tested. 
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Figure 6.- Velocity and dynamic pressure of cylinders tested. 
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Figure 7 .- Time history of fineness-ratio-2.71 cylinder. Model I 









Figure 8.- Time history of fineness-ratio-4.00 cylinder with flared 
afterbody and rounded nose. Model II. 
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Figure 8.- Concluded. 
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-1.5 -1.0 -.5 0 .5 1.0 1.5 



•5 1.0 1.5 


-1.5 -1.0 -.5 0 .5 1.0 


(a) M = 2.09 to M = 1.66. 


(b) M = 1.66 to M = I.34. 


Figure 9»- Basic data cross plots of force and moment coefficients f 
the fineness -ratio -2. 71 cylinder. Model I. The time sequence is 
indicated by the symbols O, and A. 




26 



NACA RM L58A16 




-2.5 -2.0 -1.5 -1.0 -.5 0 .5 1.0 1.5 2.0 2.5 


(c) M = 1.35 to M = 1.18. 
Figure 9- - Continued. 
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(d) M = 1.09. to M = 0.84-. 
: Figure 9.- Continued. 
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(g) M = O .58 to 


M = 0.49. 


(h) M = 0.49 to M 


0.39. 


Figure 9 .- Concluded. 
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Figure 11.- Variation of rolling velocity with Mach number 
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O Fineness -rati 0-2. 71 cylinder of this test. 
(Cc where C N < 0.1 and Cy < 0.1) 

□ Fineness -rati o-2. 56 cylinder of reference 3. 
(Cc where CijCO. 1 and Cy<0.1) 



Figure Ik.- Variation of longitudinal -force coefficient with Mach 
number . Shaded area indicates measured Cq for al 1 Cjj and 
C Y values of this test. 


NACA - Langley Field, Va. 




